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Abstract 

In this paper an analysis of the Stirling cycle in thermoeconomic terms is developed using the 
entropy generation. In the thermoeconomic optimization of an irreversible Stirling heat pump 
cycle the F function has been introduced to evaluate the optimum for the higher and lower sources 
temperature ratio in the cycle: this ratio represents the value which optimizes the cycle itself. The 
variation of the function F is proportional to the variation of the entropy generation, the maxima 
and minima of F has been evaluated in a previous paper without giving the physical foundation 
of the method. We investigate the groundwork of this approach: to study the upper and lower 
limits of F function allows to determine the cycle stability and the optimization conditions. The 
optimization consists in the best COP at the least cost. The principle of maximum variation 
for the entropy generation becomes the analytic foundation of the optimization method in the 
thermoeconomic analysis for an irreversible Stirling heat pump cycle. 
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I. INTRODUCTION 



The Stirling cycle is an important model of refrigeration systems and the recent develop- 
ments in design were proposed after the new concept of finite time thermodynamics came 
into existence . Blanchard applied the Lagrange multiplier method to find out the 

COP of an endoreversible Carnot heat pump operating at the minimum power input for a 
given heating load [4j . Several papers heve been devoted to propose mathematical functions 
to optimize thermodynamics cycles starting from different initial conditions and focusing 
on the total cost and efficiency. The definition of optimization that we adopt in this pa- 
per is the best COP at the least cost. The performance of the different heat engine and 
refrigeration systems were investigated using the concept of finite time thermodynamics, 
of the ecological approach and the thermoeconomic analysis □ D □ □ □ On the other 
hand, the key-role of entropy generation maximum has been recently demonstrated in ther- 
modynamics analysis of the irreversible processes and it has been shown that it represents 
a new criterion for determining the conditions for stability In this 

paper we show that this criterion represents the thermodynamic foundation for some recent 
results obtained in thermoeconomic analysis of the Stirling heat pump cycle. We start from 
studying the time evolution of an open system and we take as working hypothesis that it 
evolves in the optimization of the irreversibilitiy due to entropy generation. F function is well 
suited to show the system evolution related to the irreversibility due to entropy generation, 
our aim is to propose an approach based upon the natural behaviour of the thermodynam- 
ics/thermoeconomic system as a groundwork for the optimization analysis. The evolution 
of an open system is considered natural when it moves in order to get the optimization of 
the entropy generation. 



II. THE THERMODYNAMIC ANALYSIS 



The working substance of the Stirling cycle may be a gas, a magnetic material, etc., and 
for different working fluids the performance of the cycle are quite different. The Stirling cycle 
with an ideal gas consists of two isothermal and two isochoric processes. It approximates 
the expansion stroke of the real cycle by an isothermal process to whom heat is added to 
reach the temperature T c from a heat source of finite capacity whose temperature varies 
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from Tli to Ti2- The heat addition to the working fluid is thought as an isochoric process: 
heat is going towards the heat sink of finite heat capacity that gets a temperature variation 
from Tjji to Th2- The heat rejection from the working fluid to the regenerator is modelled 
as an isochoric process which completes the cycle itself. Let Q c and be the amount of 
heat absorbed from the sources at the temperature T c and T ft respectively, during the two 
isothermal processes jl[: 

Qh = C H e H {T h -T H1 )t h (1) 

Q c = C L e L (T L1 -T c )t L (2) 

where Cjj is the heat capacitance rate of the sink reservoir, Ci is the heat capacitance 
rate of the source reservoir, tn is the heat rejection time, ti is the heat addition time, 
€h is the effectiveness of the heat exchangers for the hot-side and is the effectiveness 
of the heat exchangers for the cold-side. These cycles do not possess the condition of 
perfect regeneration, hence it is assumed that the loss per cycle, AQ p, is proportional to 



the temperature difference of the two isothermal processes as follows [1 
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3: 



AQ R = nc f (l-e R )(T h -T c ) (3) 
where ct is the molar heat capacity of the working fluid and n is the number of moles. The 

n 

Gouy-Stodola theorem [20] states that the thermodynamic work burnt in the irreversibility 
due to the entropy generation is equal to the product between the lowest source temperature 
and the entropy generation, i.e. total entropy is equal to the isolated system entropy plus 
the irreversibility due to entropy generation. Considering the Gouy-Stodola theorem and 
the definition of the entropy due to irreversibility ASi rr [13j , the last one can be written as: 

AS irr = = nc/(l - e R ) ^ Th = nc f (l - e R )(x - 1) (4) 

with x = Th/T c . The theorem of maximum entropy generation states that the entropy 
generation is maximum at stationary state Q|. This theorem allows a new approach to 



irreversible processes as it is proved in a lot of different applications in hydrodynamics 1Q | , 
engineering thermodynamics jllj ]. rational thermodynamics and biophysics [3, 14 1- 
Hence applying it here, we argue that equation|3]must be a maximum in the thermodynamics 
stability: this equation described the natural behavior of the thermodynamics system. 
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III. THE THERMODYNAMICS FOUNDATION OF THE THERMOECONOMIC 
ANALYSIS 



The objective function F of the thermoeconomic optimization recently proposed is 
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Ci + C e 

with Qh = heating power, Ci and C e refer to annual investment and energy consumption 
costs, and are defined as: 

d = a(A H + A L + A R ) + b^lLZ^l (6 ) 

^cycle 

Ce = b Qh - Qc (7) 

tcycle 

where a is a constant directly proportional to the investment cost of the heat exchanger 
and is equal to the capital recovery factor multiplied by the investment cost per unit heat 
exchanger area. Ah + Al + Ar is the heat exchanger total area, with Ah the heating area, Al 
the heat source area and Ar the regenerative area, b is the capital recovery factor multiplied 
by the investment cost per unit power input and t cyc i e is defined as: 

tcycle — + + tn (8) 

with 



t R = 2a(T h - T c ) = 2aT c (x - 1) (9) 

where a is a constant that depends upon the kind of working fluid used in the cycle, 
and shows that the working time of the regenerator (a sort of recovering time towards the 
initial conditions in thermoeconomics) is proportional to the difference of temperature. In 
the thermoeconomic analysis of an irreversible Stirling heat pump cycle the function F has 
been used to evaluate the ratio of the higher and lower source temperature in order to reach 
the optimization of the cycle itself. The common solution, based upon the application of 
the variation method, consists in evaluating the maxima of F function, solving the equation 
5F = 1] applying the variational method. 
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Now, from |H1 and [7| the El becomes: 



Qh 



a{A H + A L + Aft) + {b + b' 



Qh — Qc 
i 

tcycle 



(10) 



Starting from the relations El 03 IHHTUl we can argue that the objective function F of the 
thermoeconomic optimization is related to the entropy generation as follows: 



Q 



H 



a(A H + A L + A R ) + (b + b') 



Qh — Qc 



2aT r 



t H + t L + — 

nc f (l-e R ) 

which can be easily written after few algebric operations: 

r 3 + T^ASirr 



-AS, 



(11) 



(12) 



with 



Ti= Q H (t H + t L ) 

Y _ 2aT c Q H 
2 nc f (l-e R ) 

r 3 = a(A H + A L + A R ) (t H + t L ) + (b- b') (Q h - Q c 

2aa(A H +A L +A R )T c 



4 ncf{X-e R 

From equation El we can argue that the variation of the function F is proportional to 
the variation of the entropy generation: 



with 



Then it follows that 



SF 



r 2 r 3 - r 1 r 4 
r 3 + r^ASir. 



■5(AS r , 



(13) 



r 2 r 3 ^ 



(14) 



S(AS irT ) = 0^5F = 
5 



(15) 



In this way it has been stressed the relation between the economic analysis and the 
thermodynamics. In the economical analysis the function F was introduced in several papers: 
we need to know its upper and lower limits to fulfill the basic conditions of optimization, but 
no physical explanation has been up to now given about this method. Here we prove that 
the limits of the F function are directly correlated to the entropy generation in the state 
of stability and related to the optimization of the cycle. Hence the optimization, which 
consists in the best COP related to the least cost, can be obtained in the conditions of 
natural stability for the open systems. The evolution of an open system is defined natural 
when it moves to get the optimization of entropy. The advantages of this method consist 
in exploiting the natural dynamics of the system in order to reach, following its natural 
behaviour, the optimum by the shortest way (i.e. the lower cost). 



IV. CONCLUSIONS 



The thermodynamic and thermoeconomic analysis of the optimization of an irreversible 
Stirling heat pump cycle is presented in relation with its thermodynamic foundation. We 
proved that the principle of maximum variation for the irreversible entropy is the analytic 
foundation for the optimization method recently introduced in the thermoeconomic analysis 
for an irreversible Stirling heat pump cycle. Of course it represents not only an analytical 
and mathematical groundwork, but also the physical and thermodynamic foundation for 
the method itself, as a consequence of the physical meaning of the principle of maximum 
entropy variation in thermodynamics The optimization method is a 



useful tool to design thermodynamics systems characterized by lower working costs. The 
principle of maximum variation allows a deeper thermoeconomic analysis focused on the 
stability conditions. 
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Nomenclature 



a capital recovery factor times cost 

per unit heat area 

A area [m 2 ] 

b capital recovery factor times investment cost 

per unit power input 

c molar heat capacity [Jmole^k^ 1 ] 

C heat capacitance rate \kWK~ 1 

COP Coefficient of Performance 

n number of moles [mole] 

Q heat [J] 

S entropy [JK^ 1 ] 

t time [s] 

T temperature [T] 



Greek letters 

e effectiveness 

5 differential 1 s ^l, i -§^dz i 

A finite variation 



Subscripts 

f fluid 

h sink side 

H heating 

irr irreversible which is related 

to the entropy generation 

L heat source 
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